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We analyze the Josephson current in a junction consisting of two superconductors (S) and a 
ferromagnetic layer (F) for arbitrary impurity concentration. In addition to non-magnetic impurities, 
we consider also magnetic ones and spin-orbit scattering. In the limit of weak proximity effect we 
solve the linearized Eilenberger equation and derive an analytical expression for the Josephson 
£j . critical current valid in a broad range of parameters. This expression enables us to obtain not 

only known results in the dirty and clean limits but also in a intermediate region of the impurity 
concentration, which may be very important for comparison with experimental data. 
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I. INTRODUCTION 



Negative Josephson coupling has been predicted almost three decades ago 0, 0. Bulaevskii et al [l[ obtained 
this coupling from a model of a tunnel Josephson junction containing magnetic impurities. By tunneling from one 
superconductor to the other, electrons are scattered by these impurities. Explicit calculations led the authors of 
Ref.0] to the conclusion that under certain conditions the Josephson critical current I c might change its sign. 
O \ Later on, Buzdin et al. suggested that a similar behavior might take place provided the insulating barrier was 
substituted by a ferromagnetic metallic layer (F)[U]. Since then, the 0-tt transition, i.e. the change of sign of I c , has 
^vq been the focus of study of many theoretical works (see 0, d, 0] and references therein) . 

. However, in spite of the theoretical progress, the 0-tt transition has been observed only recently in SFS junctions 
with characteristic thicknesses of the F layer of the order 100 A or less 0, [1 1, El El El El El • Although a 
qualitative explanation for the observed drops of I c as a function of thickness or temperature could be given in the 
framework of the known theories, a quantitative description is, in many cases, still lac king . This is due to the fact 
that the authors of the theoretical works considered either the pure ballistic case @, EI M EE HI or the dirty 
limit using the Usadel equation [llllll- br the pure ballistic case, the critical current I c oscillates with the period 
vp/2h as a function of the thickness 2d of the F layer, where vf is the Fermi velocity and h the exchange field in the 
ferromagnet. In the diffusive limit, both the period of oscillation and the decaying length of the function I c (2d) are 
equal to y // D/2h, where D — v f t/3 is the classical diffusion coefficient and r is the momentum relaxation time due 
p ■ to the impurity scattering. 

The diffusive limit is realized provided either the exchange energy h is not too large or the mean free path I is short 
enough. More precisely, the condition hr << 1 should be fulfilled. However, in many experiments this is not the case. 
In particular, in those performed with strong ferromagnets d El, E3 • Theory considering an arbitrary value of hr 
£j . & n d / has been presented in Ref. (2lj , where the Eilenberger equation for the condensate has been solved in the weak 
proximity effect limit. In particular, it was shown that in a quasiballistic case (or in case of strong ferromagnets), i.e., 
. Sh ' when the condition 
K> ■ 

' hr»l (1) 

c3 ; 

is fulfilled, the critical current oscillates with the period vp/2/i and decreases exponentially over the mean free path 

I = VpT. 

In all theoretical works mentioned above the ferromagnet was modeled as a normal metal with an exchange field 
acting on the spin of the conduction electrons, while depairing factors such as spin-dependent (SD) scattering on 
magnetic impurities, condensate flow (due to an internal magnetic field in F) and spin-orbit (SO) scattering have not 
been taken into account, although they may play an important role in the proximity effect. 

The effect of the spin-orbit scattering on the critical Josephson current in SFS junctions in the diffusive limit has 
been studied in Ref. [22|, [23| and in a recent work [24], where it was shown that the SO interaction affects the 
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characteristic length of the decay of the critical current with increasing the thickness of F. Also in the diffusive limit 
the effect of the scattering on magnetic impurities was considered in Refs. [24l. [25j. 

Less attention has been paid though to the study of depairing effects in the quasiballistic regime when the condition 
(TJJ) is fulfilled. An attempt to solve the Eilenberger equation with account for the SO and SD scattering and for 
arbitrary I was undertaken in Ref.[26|. However, the solution suggested by the authors of the latter work is valid 
only in the diffusive limit. As we will see, in the limit determined by Eq.([T]), the expression obtained in Ref. [26| 
is a small part of the total solution. The SD scattering was also considered in a recent work (2?| in the case of a 
weak proximity effect and arbitrary mean free path Unlike Ref. (2lj . the authors suggested an approximate solution 
for the Eilenberger equation and neglected the SO scattering. Therefore, the problem of calculating the Josephson 
current Ij in a general case of arbitrary mean free path I taking into account different depairing mechanisms remained 
unsolved. 

In this paper, we attack this problem calculating the current Ij through an SFS structure in the general case of 
an arbitrary mean free path I taking into account diverse depairing mechanisms. Following the method presented 
in Ref.[2l| we solve in the next section the Eilenberger equation for the case of the spin dependent scattering on 
magnetic impurities. As in Refs. [U [H, [13]) we assume that the proximity effect is weak, i.e., that the amplitude of 
the condensate function induced in the ferromagnet is small. This assumption is reasonable even for low temperatures 
due to the large mismatch of the electronic parameters of the superconductor and ferromagnet resulting in a strong 
reflection at the S/F interfaces. Using the exact solution obtained for the condensate function we derive a general 
expression for the Josephson critical current through the SFS system. This expression can be used for calculating 
the current at arbitrary impurity concentration. In particular, we calculate I c in the diffusive and the quasi-ballistic 
limit. In section III we consider the effect of the SO scattering on the critical current I c for these two cases. 



II. SOLUTION FOR THE EILENBERGER EQUATION AND DERIVATION OF THE JOSEPHSON 

CURRENT 

We consider an SFS layered structure. The thickness of the F layer is 2d and the F/S interfaces are located at 
x = ±d. The thickness of the S layers is assumed to be infinite. Scattering of electrons by magnetic impurities in a bulk 
superconductor has first been studied by Abrikosov and Gor'kov [28| using microscopic Green's functions. For non- 
homogeneous finite systems it is more convenient to use quasiclassical Green's functions determined by the equations 
derived by Eilenberger, Larkin and Ovchinnikov [29l . [30| . In order to justify the applicability of the Eilenberger 
equation, we also assume that the distance between the superconductors, i.e. the thickness of the F layer, is larger 
than the mean free path I (see Refs. [3l|, |32().We consider the case of a weak proximity effect, i. e. when the amplitude 
of the elements of the condensate matrix function / is assumed to be small: |/| << 1. In this case one can linearize 
the Eilenberger equation that in the presence of the scattering on non-magnetic and magnetic impurities takes the 

form @, i m, m 

sgn W f 3 ZeV/± + (nf)± = (1 - 2X Z - A ± )(/±) + A x </=,=) (2) 

where T3 is the z-component of the Pauli matrices f, e = vf/vf is unit vector in the direction of the Fermi velocity, 
k± = 1 + 2(\lo\ ± ih u )rt, uj = uj n = nT(2n + 1) is the Matsubara frequency h w = fegnw, rf 1 = t^ 1 + is the 
total scattering rate and tjv is the momentum relaxation time due to scattering by nonmagnetic impurities. The rate 
t m\ = T M / (1 + a ^ + a ±) is the total scattering rate due to scattering by magnetic impurities. The parameters a z , a± 
characterize the spin-dependent scattering. The potential of interaction with magnetic impurities can be written in 
the form 

U(r) = U (r) + U s (r)S-a/S, 

where the potential U (r) describes interaction with a spin- independent part of magnetic impurities. The impurity 
spin S is assumed to be classical. The coefficients a z and a± are given by a z — [\Us\ 2 /\U \ 2 }(S Z } / S 2 and a± = 
[\Us\ 2 l\U \ 2 ]{S 2 ) / S 2 + [\Us\ 2 /\U \ 2 ](S 2 ) / S 2 . These coefficients are related to a spin-dependent scattering rate t" 1 
used in Ref. [24 ] ■ For example, a z = 2(tm / T m ) (S 2 ) / S 2 . The angle brackets denote averaging over angles. The 
coefficients A z ,Aj^ are defined as: X z = {ct z /TMt)/{i'^j t +T^ N 1 ) and Aj_ = { a ±/ T Mt)/{T M \ +t^ 1 ). Note that a z ± << I 
& 

In the general case the condensate function / is a 4 x 4 matrix in the particle-hole (Gor'kov-Nambu)£g)spin space. 
However, in the case of a homogeneous magnetization considered here and in the absence of the spin-orbit scattering 
(the SO scattering will be taken into account in the next section), the function / is diagonal in the spin space. Then, 
the function /+ in Eq. ([2]) is defined as / + = f aa with a = 1 (a is the spin index), while the other diagonal element 
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is given by /_ = $22 = ~f+(~h)- I n the present case of a planar geometry, the function / depends only on the 
coordinate x. 

It is convenient to represent / as a sum of a symmetric and antisymmetric part with respect to the momentum 
direction: f±(x) — s±(x) + a±(x). As follows from Eq.Q, the antisymmetric part d{x) is related to the symmetric 
one by the expression 



a± — — sgnw (pd/n±) T^ds±/dx 
while the symmetric part s(x) in the ferromagnetic region (— d < x < d) obeys the equation 

fi 2 l 2 d 2 s±/dx 2 - k 2 ± s± = -K±(h) [(1 - 2X Z - X±)(s±) + X±(s T )} 



(3) 



(4) 



where I = vpr t is the mean free path and /i = v x jvp — cos9. These two equations should be complemented by the 
boundary condition [3~ij 



--±d= T7(/^)sgnw (t 3 /s) ,fi>0 



(5) 



where j(^) — T(/x)/4 and T(/i) is the transmission coefficient. The latter is assumed to be small and therefore the 
matrix f s in the superconductor has its bulk form f s (±d) = f s (f2 cos(<p/2) ± fi sm(<p/2)), where f s = A/y/oj 2 + A 2 . 
Note that the boundary condition Eq. does not take into account spin-flip processes at the interface. Boundary 
conditions for magnetically active interfaces were derived in Ref. (35j and used in Ref. 36] by calculating the 
supercurrent through a pure ballistic half-metallic layer. 

In order to derive an expression for the Josephson critical current avoiding straightforward but cumbersome cal- 
culations we first neglect the terms proportional to A j_ . This term will be included again at the end of this section 
where we will present the final expression for the current. The way how to obtain the solution of Equations (J^HHl) wa s 
presented in Ref. [2lj . One formally extends the solution over the whole x-axis, performs a Fourier transformation 
and obtains finally the transformed function §k± from an algebraic equation. Following this scheme we obtain an 
expression describing the spatial dependence of the condensate function 



s±(x) = ± 



dk 2k±1 
2^M fe± (/i) 



E 



f St n exp(ikd(2n + 1) — ikx) 



K±(l 



2A z )iV ± 1 



M7(m) 
A4±(/i) 



(6) 



with f s 



/ s [f 2 cos(<^/2) + {-l) n n sin(¥>/2)], M k± (fi) = (M) 2 + 4 and N± = l-(1- 2A 2 ) K± (Af fe± 1 (M)> 



The average n(M k can easily be calculated (we drop the subindex ±) 



kl 



■ tan 



kl 



1 in 

TT7 ln( - — 
likl ik 



kl, 



kl' 



(7) 



The x-dependence of the condensate function, in particular, its exponential decay is determined by the singular points 
of the integrand in the complex k plane of Eq. From Eq. ([7]) one can immediately see, that the function in the 
integrand has poles at kl = ±m/|/x| and branch points at kl = ±ik±. The poles lead to the exponential decay of s(x) 
like s po i(x) ~ Cpoi(fi) exp(— n(d + x)/l\/j,\) (in the vicinity of the left superconductor), whereas the branch points lead 
to terms in the solution that decay as Sbr(x) ~ Cbr(fJ-) exp(— n(d + x)/l) (where —d < x). This means that in the limit 
Tt < 1 the terms s po i (x) exponentially decrease over an angle-dependent distance of the order of the mean free path 
/. On the other hand, the terms Sb r (x) decay exponentially over an angle-independent distance of the order I. Both 
the terms, s po i(x) and Sf, r (x), oscillate with the periods h/vp\fi\ and h/vp, respectively. The amplitudes C po ;(/x) and 
Cb r (/x) depend on the parameters of the system. 

In the general case their calculation is rather complicated and a numerical analysis is needed. However, as we 
will see in the next paragraphs, in the quasi-ballistic and diffusive cases one can obtain explicit expressions for the 
condensate function s(x). If one is interested in an intermediate case it is convenient for numerical calculations to 
perform the integration over momentum in Eq. ([B]) taking into account the relation (see, for example [37| ) 



exp(i2nfcd) = — S(k — k m ), (8) 

n— — 00 

where k m = irm/d. 



d 

m— — 00 
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Writing the condensate matrix as s± = sfr\ + sff2 one finally obtains 



± ,1> K ± exp[i(x/d — l)(2m + 1)tt/2] 



— 



m— — oo 



E 



exp[i(x/d — l)m7r] 



m— — oo 



Af. 



K ± {l-2\ z ){NfJ 



K± (l-2A z )(iV 2 ± m+1 )- 1 



M7(m) 

M 2m+lO) 



M7 



M7 



/ s cos y>/2 



/ s sin<^/2 (9) 
(10) 



where = (filnm/2d) 2 + k± and JV^ = 1 — (1 — 2X z )k±(1/M^) . Knowing the condensate function induced in the 
ferromagnetic region one can calculate the Josephson dc current Ij through the SFS junction. This current is given 
by 



I = {iriT/4)(e 2 /h)(k 2 F S/ir 2 )Tr 



t 3 ®cto E (M7[s(rf),/ S (<i)]-) 



(11) 



where [s, / s ]_ = sf s — f s s, ctq is the unit matrix and the symbol Tr stands for the trace over the 4x4 matrices (see 
Using Eq. ([6]) , one can write the current in the form I = I c sin ip with 



°° / + + / ' ( i '1 

J c = A(2*T)Re £ /'(-i)^^) ' - 2A.) ' '' (/ ' ' 

a Mm {[I) l\ m 



m— — oo,uj>0 



' M7(/")] 



(12) 



where A = ((e 2 /K)(k 2 F S/ir 2 ) and M m (fj.) = M+(/x). Eq. (fT2")) is the most general expression for the Josephson current 
in terms of the solution of the Eilenberger equation. In the next sections we give expressions for the critical current 
I c in the quasiballistic and diffusive limits. 



The quasiballistic case: \k\ » 1 

The quasiballistic case corresponds either to a strong ferromagnet (h > r _1 ) or to a clean sample (Tt > 1). As 
follows from Eq.©, in this case N « 1 and the second term in the square brackets in Eq.© is much larger than the 
first one. Calculating the residue at the pole kl = in/\n\, we obtain for s{x) = s + (x) 



s(x) = 7(/i)/ s 



A cosh( x / L a b(n)) , smh(x/L a h(u)) . , 

f2 ■ uV// cos ^/2 + h .),' qb Y\ sin y/2 
smh(d/ LqbifJ,)) cosh(d/L qb ((j,)) 



(13) 



with the length L q b(p) characterizing the quasi-ballistic case defined as L q b(/i) = l\fi\/n. 

It is seen for example that in the vicinity of the left superconductor, i.e., at 1 << (x + d)/l << d/l the condensate 
function s(x) oscillates with the period nvp/h\fi\ decaying over the angle-dependent mean free path l\fj,\ : s(x) ~ 
exp(— kx/l\fi\) 21]. Thus, if the exchange energy is large, h > r _1 , the function s(x) experiences many oscillations 
over the decay length I. Note that the period of oscillations and the decay length strongly depend on the angle 
H = cos 9. 

This result contradicts the conclusion of Ref. [26| where the solution for Eq. ^ was taken in a form of an exponential 
function with an angle- independent exponent. As can be understood from Eq.©, the solution obtained in Ref. [26j 
is not general and corresponds only to a contribution from the branch points, i.e., from the first term in the square 
brackets in Eq.Q. However, in the limit of large k this term is small compared to the other one (see [HI). 

We see from Eq. ((4]) that in the quasiballistic case the depairing leads only to a renormalization of the scattering 
time t — > r t . 

The main contribution to the critical Josephson current in the quasiballistic regime stems from the second term in 
the square brackets in Eq. (fT2"j) . Calculating the residues at the poles of k 2 , (//), we obtain (cf. Ref. [2l[) 

*-^£^(U%Sw»), (14) 

If the thickness of the F layer considerably exceeds the value of L q b(l), one obtains for the critical current I c 

J, = A(2nT) jr /Ml) *^/ v VF) exp[-(2d/0(l + 2cr t )] (15) 
uj>o ' Vf 

Eqs. (|14j|15|) fit well recent experimental data [1, [l2j 
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The diffusive case: \k\ ~ 1. 



In the diffusive limit the conditions hr, Tt << 1 hold and the characteristic length of the spatial variation of s(x) 
is much larger than the mean free path I. Therefore, we obtain from the following expressions for the parameters: 



x dif 



(1 - 2X z )(kl) 2 /3 and k 2 « M k (p) « 1 (c/ Eq. ©). 



The behavior of s{x) at distances from the superconductors larger than the mean free path I is determined by the 
residue of the pole in Nk- Equation ^ finally yields 



s(x) 



y/3(l - A z ) 

Kdlf 



. . . \ cos(y>/2) 
smh^d/Ldi/j 



sinh(a;/L dl /) . 



cosh(d/L 



dif) 



(16) 



where k^* — X z + 2(\u\ — ih u )T t , and the characteristic length Ldif of the condensate decay in F is Ldif — 

1^(1 - 2X z )/3K 2 dif « l/(K dif VS). 

If the exchange energy h is much smaller than the spin-dependent scattering rate, hr t << X z , we obtain for the 
characteristic length Ldif ~ l/y/3X z = \J Dr sp /2, where r sp = t/X z . Thus, as expected, the decay length in this 
case is related to the spin-dependent relaxation time (see, for example , [391) . In the opposite limit, h >> l/r sp , 
the characteristic length Ldif is given by the well known expression [H. Il9l. l2~il [25| : Ldif ~ — i)\fihrt\ = 

(1/2) (1 + i)y/D/h. Note that in this case the condensate function both decays and oscillates on the same length 
jLTjh. 

If one takes into account an internal magnetic field B inside the ferromagnet given by B 
in the limit hr t 



« X z is equal to: Ldif = (2dB/4> ) 
6q = irhc/e is the magnetic flux quantum (see, e.g., [40(). 



(Dt sp /2Y 



4tvM, the length L d if 
where M is the magnetization in F and 



In the diffusive limit, the critical Josephson current I c is determined by the first term in the square brackets of Eq. 



(Q2J), i.e., by the poles of the function AT, . Then, we find 



A(2irT) V3/ s 2 Re 

uj>o 



K dif smh(2d/ Ldif) 



(17) 



where Ldif is defined in Eq.p^]). 

An alternative way to obtain Eq. <\17\i is to solve the Usadel equation, which was done in many publications 0, fl9l . 

MM- 

For an arbitrary value of the parameter hr t the critical current can be computed from Eq. (|12[) . In Fig. 1 we show 
the dependence of the absolute value of I c on the thickness 2d of the ferromagnet for different values of hr t taking 
into account only the scattering by non-magnetic impurities, i.e when — > 0. One can readily observe the crossover 
from the diffusive case (hrt — 0.2), when the decay length and the period of the oscillations are the same WD Jh, to 
the quasiballistic case (hr t — 2.2), when the period of the oscillations is irvp/h while the decay length is of the order 
of the mean free path I. 

To simplify numerical calculations, we assume that the transmission parameter 7(/i) is peaked at y, = 1 and replace 
it in Eq. (|12|) by a delta-function, i.e. we assume that only electrons with momentum direction perpendicular to the 



S/F interface are transmitted (the angle-dependence of "the transmission coefficient" 7(/x) depends on properties of 
the S/F interface). Another limit ('y(fj) = const) was assumed in Refs.[9l. [l2j where the expression for I c derived in 
Ref. [2l| in the absence of depairing mechanisms was used for comparison between theory and experimental data. In 
Fig. 2 we represent the critical current as a function of 2d for three different values of the parameter hrt- One can 
see that in all cases an increase of X z leads to a decrease of the amplitude of the condensate in the ferromagnet. 

Finally we write down the expression for the critical current for the case when perpendicular fluctuations of the 
exchange field are taken into account, i.e., when is finite. In that case we obtain for the critical current a rather 
cumbersome formula 



I c = A(2nT)Re £ / S 2 (-1)"( M7 ( M ) 



2k+ 



n= — og.lj >0 



(1 -2A- - A ) j h + X-(J^^ 



d M+{u) (\N+\ 2 - X\\k + \ 2 |((M+)-i}| 2 ) ' 

- a ± | K+ i 2 <(M+r v)) M , rv)>, 



-A i K-N, 



A ± | K+ | 2 <(M")- V)L (»W)(M+)-\n')). 



^M7(M)>, (18) 
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hx =2.2 



hx=1.2 



hx=0.2 




2d/l 



FIG. 1: Dependence of the absolute value of the Josephson critical current I c on the thickness 2d of the ferromagnetic layer for 
different values of hr. Here At = 0.1 and Tr = 0.05. Only scattering by non-magnetic impurities is considered. 

This equation can be evaluated only numerically. We represent the function I c (d) for different values of Aj_ in Fig. [3] 
Note that, again, with increasing Aj_ the amplitude of the condensate decreases. This is clear from the physical point 
of view because any depairing factors lead to a suppression of the condensate amplitude. 



III. SPIN-ORBIT SCATTERING 



In this section we consider influence of the spin-orbit (SO) scattering on the Josephson current. For simplicity we 
neglect the spin-dependent scattering analyzed in the preceding Section. In the presence of the SO scattering, the 
condensate function is no longer diagonal in the spin space. Therefore, instead of the 2x2 matrix /, we have to 
introduce a more complicated 4x4 matrix / (see for example @). 

The Josephson current in SFS junctions in the presence of the SO interaction was analyzed in the diffusive limit in 
Ref. [22|. Here we focus on the opposite quasiballistic limit. 

The linearized Eilcnbcrger equation for / has the form 



sgnujT 3 leWf + kf = p(f) - A so (/) so , 



(19) 



--1 -1 | (O /QW-l 

vviicnj ; aiiu i so a 

scattering and the spin orbit scattering time, p = T t /r w 1, \ so = T t /r so << 1. 
ihuU^jT-t, which means that k + in Eq.([2]) is the (k)h element of the matrix k. 
averaging: 



In Eq. (|19[) . I = vpT t ,r t 1 = r 1 + (2/3)t so 1 , where r and t so are the momentum relaxation time due to a potential 

The matrix k = n(h) = 1 + 2(|w| + 
The angle brackets mean the angle 



</) = (1/47T) / dSlfiQ), (f) so = (1/47T) 



(S x e)i/(fi')(3 x e ) A 



(20) 



where the vector S has components Si: S = (<7i, <T2, 03 <g> fa). 

As before, we represent / as a sum of the antisymmetric and symmetric parts: / = a + s, where the antisymmetric 
part a is expressed in terms of the symmetric function as a — —k^ 1 sgnajT 3 eilds/dx and the symmetric part obeys 
the equation 



eil'd's/x 2 + k 2 s = k[p(f) - A so (/) so + 2(eiZ) 7 



f Sin S(x - d(2n + 1))] 



(21) 






FIG. 2: Dependence of the absolute value of Josephson critical current I c on the thickness 2d of the ferromagnetic layer 
different values of hr and \ z . Here Art = 0.1, Tr t = 0.05 and X x = 0. Without loss of generality we have set r^ 1 — > 0. 
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FIG. 3: Dependence of the the absolute value of Josephson critical current I c on the thickness 2d of the ferromagnetic layer 
for hr = 1.2 and different values of A_i_. Other parameters are taken as in Fig. [1] 



Taking into account the structure of the f s matrix in the spin space we see that the matrix / s>n coincides with the 

one presented above, f s<n = u 3 f s ^ n . 

Our task now is to solve Eq. (j2~lj) . The presence of the term of the SO scattering makes this task more difficult 

than previously. In order to simplify the problem, we use the usual smallness of X so [281 ]. An additional simplification 

comes from using the quasiballistic case when the value of |k| is large. 

In order to find the solution of Eq. (|21[) we represent the Fourier transform s k in a form of an expansion in the 

small parameter A so : s k = S k + S Sk- 
in zero order approximation in A so and in the main approximation in the parameter \k\ we obtain 

oo 

S k = 2kle ll M-\e 1 ) /»,n exp{ikd(2n + 1)) , (22) 

n— — oo 

where M k (e{) is a 2x2-matrix in spin-space with component (1,1) equals to M k+ {e\) and component (2,2) to M k +(ei). 
Although the matrix Sk is diagonal in the spin space, the correction 5S k is not. It is equal to 



SS k = -kM^ityXsoiSkifl')}^ (23) 
Using Eas.([2"2"j). (12TJ1) , one can represent the average {S k (0,')) so as 

(5fc(n')).o = (Sko(tt') (8 a A - S k3 (il') ® a 3 A 3 + 2f 3 S k0 (n') ® HfrAj. - a 2 A 2 )) (24) 

where A Q = A (Cl,n>) = e f (e§-el)+e?(e|-e?)+ef (e?+ef); A 3 = A 3 (n,Cl>) = ef(el+e|)+e?(e|+e?)+ef (e?+e|); 
A\ = e±e2e 3 ; A^ = e 2 eie 3 . The matrices S k o, 3 are defined with the help of the relation Sk = 0q <S> S k o + 03 ® S/C3- 
Using Eq.([21]) we obtain 

Sko = 2lex^iF s Im-^-;S k3 = 2l eil F s Re-^- (25) 
M k+ M k + 

where F s = YlnL-aa exp(ifcd(2n + 1)) and the matrix / Sin has been introduced in Eq.©. 

Note that due to the last two terms in Eq. ([24|) the condensate matrix SS k contains not only the triplet component 
with the zero projection on the z-axis but also the triplet components of the type t1\ 11- However, in the lowest order 
in A so these components do not contribute to the Josephson current because they are odd functions with respect to 
the inversion e± t 2 =>■ — ei ( 2- 
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The correction 5I C to the Josephson current due to spin orbit scattering is given by the expression 

dk 



SI C = A(2niT)Tr 



h <S> <5o ^2 

LU— — QG 



2tt 



(e 11 [SS k (n)J s (d)]-) 



exp(-ikd), 



Finally, we obtain the correction SI C to the critical current originating from the SO scattering 



I 



5I c = A(2kT)- fU{eie[ 1 (n) 1 (n')[-A Im--t-Im 



M 2 + n (0) M+(fi') 



A 3 Rc 



-Re- 



M+(0) MU^Y /a 



(26) 



(27) 



The structure of this equation is similar to the one of the contribution of the first term in the square brackets in 
Eq. (fT2f to the Josephson current. This contribution is also small in comparison with the current I c determined by 
Eq. (fT4")) with L q b(n) — 2|mI/ k arL d ' = v ft~u where t^ 1 = r _1 + (2/3)t so 1 . This small correction to the critical current 
I c can be calculated numerically. In the diffusive case the SO scattering has been studied in Refs . [a. 1X91. \24. |25|| . 



IV. CONCLUSIONS 



Assuming a weak proximity effect we have derived the exact expression for the Josephson current through an SFS 
junction for arbitrary impurity concentration and in the presence of spin-dependent scattering. In the quasiballistic 
and diffusive limits this expression takes a simple form. In the former case, the parameter l/(/ir) is small. In 
the main approximation the expression for the critical Josephson current is reduced to Eq. (|14[) that agrees with 
the corresponding equation of Ref.(2l| provided the momentum relaxation time r is replaced as: r _1 =>• r t _1 = 
(1 + A^)r _1 + (2/3)r s ~ 1 . Therefore, the deparing leads only to a renormalization of the mean free path I that 
determines the decay of the condensate function in F and of the Josephson critical current I c . Oscillations of these 
quantities have the period irvp/h. 

In the diffusive case the oscillation period and the decay of the critical current I c are determined by the value of 
the product tiTd ep where the time Td ep is defined as Tdep = min{r sp , t so }. If the exchange energy lies in the interval 
T dep < h < t -1 , then the period of the oscillations of I c is 2-n^J Drd ep / (hTdep) and the decay length is ^ Drdep ■ lu 
the limit h < the period of oscillations and the decay length are determined by the value yjDjh. The general 
expression, Eq. (112j) . may serve for numerical calculation of the critical current in the intermediate region of parameters. 
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